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(N . ABSTRACT 
X) \ 

^ . The apparent absence of stars in the Milky Way halo with —5 < [Fe/H] < —4 

^ . suggests that the gas out of which the halo stars were born experienced a period of 

^ . low or delayed star formation after the local universe was lit up by the first, metal- 

^ ! free generation of stars (Pop III). Negative feedback owed to the Pop III stars 

! could initially have prevented the pre-Galactic halo from cooling, which thereby 

0> ■ delayed the collapse and inhibited further star formation. During this period, 

^ ! however, the nucleosynthesis products of the first supernovae (SNe) had time to 

O ! mix with the halo gas. As a result, the initially primordial gas was already weakly 

Q ! enriched in heavy elements, in particular iron, at the time of formation of the 

r-| ! Galactic halo. The very high, observed C/Fe ratios in the two recently discovered 
hyper metal-poor stars ([Fe/H] < -5) HE 0107-5240 and HE 1327-2326 as well 

O \ as the diversity of C/Fe ratios in the population of extremely metal-poor stars 

\ ([Fe/H] < —3) are then naturally explained by a combination of pre-enrichment 

\ by Pop HI stars and local enrichment by subsequent generations of massive, 

.!-H \ rotating stars, for which the most massive ones end their lives as black hole- 

^ \ forming SNe, only ejecting their outer (carbon-rich) layers. The possible existence 

5^ \ of populations of mega metal-poor/iron- free stars ([Fe/H] < —6) is also discussed. 

Subject headings: Stars: abundances - Stars: individual (HE 0107—5240, HE 
1327—2326) - Stars: Population II - Galaxy: evolution - Galaxy: formation 



1. Introduction 

The unique chemical abundance patterns imprinted in the old, extremely metal-poor ([Fe/H]^ < 
—3, Beers & Christlieb 2005), stellar population of the Milky Way's halo play a key role in 



""^[A/B] = log(riA/".B) — log(?^A/"■B)0J where nj^ and n-g are the number densities of element A and B, 
respectively, and where log(nA/'iB)0 denotes the logarithmic abundance ratio of A and B in the Sun. 
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understanding the formation of the Galaxy and the nature of the first stars. During the 
formation epoch of this extremely metal-poor population, the presence of carbon, as well as 
oxygen, was crucial for low-mass star formation to occur as these elements act as efficient 
cooling agents in interstellar gas (Bromm & Loeb 2003). A curiously large star-to-star scat- 
ter is also found in, e.g., C/Fe for the stars in this metallicity regime. Two stars stick out in 
particular; HE 0107-5240 (Christlieb et al. 2002) and HE 1327-2326 (Frebel et al. 2005). 
They show the hitherto highest reported enhancements of carbon of [C/Fe] ~ 4 and are 
both depleted in iron a factor of about 20 with respect to the next most iron-poor stars at 
[Fe/H] = —4. As such, they are defined as carbon-enhanced ([C/Fe] > 1), hyper metal-poor 
stars ([Fe/H] < -5, Beers & Christheb 2005). 

Several scenarios trying to explain the individual chemical abundance patterns of HE 0107—5240 
and HE 1327—2326 have been discussed in the hterature including accretion of interstel- 
lar matter, accretion of processed material from a binary companion, self-enrichment, and 
prc-enrichment of the molecular clouds out of which the stars were born (see, e.g.. Beers & 
Christlieb 2005 and references therein). In particular, the black hole-forming supernova (SN) 
models by Umcda & Nomoto (2003) can successfully fit most of the relative abundance ratios 
in both HE 0107-5240 and HE 1327-2326 (Iwamoto et al. 2005). This, however, requires 
fine-tuning of the parameters controlling the amount of mixing and fall-back in the models. 
More importantly, although the metal-poor tail of the observed Galactic halo metaUicity 
distribution suffers from small number statistics, the presence of stars below [Fe/H] = —5, 
or rather, the apparent lack of stars in the range —5 < [Fe/H] < —4 is intriguing and cannot 
be explained by the models of Umeda & Nomoto (2003). 

Below, we will explore an alternative formation scenario for the hyper metal-poor stars 
and argue that star formation was suppressed during the initial phase of the collapse of the 
Galactic halo. Such a suppression would, as we shall see, lead to a gap in the halo metaUicity 
distribution and retain a large scatter in C/Fe, even for unevolved, single, extremely metal- 
poor stars. 

2. Star Formation and Chemical Evolution in the Early Galaxy 

Star formation in the pre-Galactic halo could have been inhibited by negative feedback 
mechanisms by the first (Pop HI) stars. Lyman- Werner radiation (dissociating H2), ultra- 
violet radiation (ionizing H), and Pop HI SN explosions (heating and removing gas) could 
initially have prevented the halo from cooling, collapsing, and forming new stars (see e.g.. 
Reed et al. 2005 for a discussion on feedback effects). Such an early feedback phase has 
also been proposed as a solution to the so called angular momentum problem in simulations 
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of disk galaxy formation (e.g., Sommer-Larsen, Gotz, & Portinari 2003). Hence, when the 
collapse of larger structures eventually occurred and the gas became re-accessible to star 
formation, the Galactic halo had already been enriched in metals, e.g., iron, by the first 
generation of SNc. Stars formed in the beginning of the second star formation epoch may 
therefore have contained small traces of iron, to a level of [Fe/H] < —5. Concurrently, ejecta 
from the second generation of SNe were dispersed in the interstellar medium (ISM) of the 
halo. However, due to the increased star formation rate (SFR), this material did not have 
time to mix as much with its surroundings before locked up in subsequent generations of 
stars. Later stellar generations should therefore display significantly higher abundances of 
iron, i.e., [Fe/H] > —4. We shall now quantify this idea. 

The stochastic chemical enrichment model used in the present study is described else- 
where (Karlsson 2005; Karlsson & Gustafsson 2005). Here we will only briefly mention 
the main features of the model including new improvements. Since observational data on 
oxygen are uncertain (Garcia-perez et al. 2006) and generally less available for stars below 
[Fe/H] = —3, we will, apart from iron, focus on the evolution of carbon. 

The ejecta of each exploding star are assumed to be evenly spread within a continu- 
ously growing volume T4iix(^), the mixing volume. The mixing volume growth is caused, e.g., 
by turbulent motions and cloud-cloud collisions in the ISM and is modelled as a random- 
walk process. Hence, Vmixit) = 47r x (o'mixt)^^^/3, where the diffusion coefficient is set to 
<7mix = 1-12 X 10"'^ kpc^ Myr~^, within the limits estimated from data given in Bateman & 
Larson (1993). Star formation is assumed to be spatially uncorrclatcd, i.e., ip = ip{t), im- 
plying that exploding stars are randomly distributed in space. With this assumption, the 
probability of finding a region enriched by the ejecta of k stars is given by the Poisson dis- 
tribution wisM{k,t) — e~^(*)//(t)'^//c!, where //(t) is the average number of stars contributing 
to the enrichment in a random point in space at time t. 

In order to properly model the build-up of carbon in the ISM, the contribution from 
intermediate mass stars should be included, taking into account the delay due to non-zero 
stellar lifetimes. Hence, the SN rate i^sn should, in all relevant equations, be replaced by the 
"stellar enrichment rate" user, defined as the total stellar death rate at time t. The death 
rate of stars of mass m, Um, is given by 



where 0(m) oc m (Kroupa 2001) is the initial mass function (IMF), iplt) is the SFR 
(Figure 1), and grim) is the lifetime of a star of mass m. Hence, 




(1) 
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mser(^) = j Um{t)dm, (2) 

9r\t) 

for t > griiTiu) and zero otherwise. The lower mass hmit is determined by the least massive 
star able to enrich the ISM at time t while the upper mass limit = 60 is fixed. 

Relaxing the instantaneous recycling approximation we need to calculate a mixing mass 
distribution (for details, see Karlsson 2005) for each stellar mass m, i.e., 



fM,M = 4/ a^{t)V^Urv)wiSM{k -l,t) 
Jadm 

X Um{t - Tv)^{t)dtdTv, (3) 

where the total supernova rate i^sn is replaced by Um- The function ai,{t) denotes the 
fraction of still-surviving stars formed at time t and cj^ is a normalization factor. The above 
probability density function is a generalization of the mixing mass distribution Jm,, 

derived in Karlsson (2005) and simultaneously gives the stellar masses (i.e. the yields) and 
the mixing masses (i.e., the amount of dilution of the ejecta) of k enriching stars. 

As an efficient cooling agent, carbon presumably played an active role in the formation of 
extremely metal-poor stars by promoting fragmentation of interstellar gas clouds. This was 
accounted for by suppressing low-mass star formation in gas with [C/H] < —3.5 (Bromm & 
Loeb 2003). Furthermore, since SN-induced star formation partly occurs in the ISM, this was 
introduced by increasing the the probability of forming low-mass stars in close connection to 
SN remnants. A minimum mass of M^ix = 3 x 10^ A4.Q was, however, assumed to be swept 
up by the remnants before stars were allowed to form (Ryan, Norris, & Beers 1996). The 
inclusion of SN-induced star formation is not crucial for the conclusions presented here. 

Carbon yields by Meynet & Maeder (2002) derived from rotating intermediate- and 
high-mass star models (metaUicity Z — 10~^) were used. Iron yields for SNe in the mass 
range 13 < tii/Mq < 30 {Z = 0) were taken from Umeda & Nomoto (2002). SNe with 
extensive fall-back, i.e., black hole-forming SNe, were assumed to form in the mass range 
30 < m/M-Q < 60 and their iron yields were set to zero. The small amount of iron supposedly 
produced in 0-Ne-Mg core-collapse SNe (8 < m/ J^q < 10) was, however, taken into account 
(Mayle & Wilson 1988). Yields of pair- instability SNe were not included (however, see § 4). 

The effect of the primordial stellar feedback was realized by assuming an exponential 
decay of the local (Pop III) star burst to a low value (cf. Sommer-Larsen et al. 2003), followed 
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by an increase of the SFR to a value of 0.4 kpc""^ Myr~^, appropriate for the Galactic halo 
during this period of time (Samland, Hensler, & Theis 1997). The star burst was normalized 
to 0.12 SNe per kpc^ over 100 million years (25 > z > 17). This particular value is arbitrary 
but lies in between different estimates of the early cosmic SFR (e.g., Springel & Hernquist 
2003; Daigne et al. 2004). The adopted SFR is shown in Figure 1. The dashed line at 
z — 17 (±4, 68% confidence errors) marks the redshift of reionization as suggested by the 
WMAP data (Spergel et al. 2003) while the line at 2; = 5 marks the adopted redshift of 
formation of the Milky Way galaxy. These redshifts were used to estimate, respectively, the 
end of the first star formation epoch (Pop III) and the beginning of the second, main star 
formation epoch in the Galactic halo. As we will see below, an initial star burst followed by 
a quiescent period of low star formation will generate a " pre-enrichment" -peak in the halo 
metallicity distribution. Such a peak would be absent for non-bimodal (e.g., constant or 
increasing/decreasing) SFRs (see Karlsson 2005, Figure 6). The particle density of the ISM 
was set to 0.1 cm~^. 



3. Results 

Figure 2 shows the predicted halo metallicity distribution. Owed to the bimodality of the 
SFR, a small group of hyper metal-poor stars appears in the range —5.5 < [Fe/H] < —5 
(implying a mixing mass of Minix ~ 10^ Mq). These stars were formed out of pre-enriched 
gas in the beginning of the second star formation epoch ai z = 5 (cf. Umeda & Nomoto 2003). 
Depending on whether the ISM at that point was, on average, enriched by 1 SN or many SNe, 
stars in this group would either be found in a few narrow peaks (the present model) or in 
one broader peak. The predicted number of stars at [Fe/H] < —5 is in good agreement with 
observations only if suppression of low-mass star formation below [C/H] = —3.5 is taken into 
account (Bromm & Loeb 2003). This suggests that the primordial/early IMF differed from 
that of the present-day, otherwise hyper metal-poor stars would have been more abundant. 

The predicted distribution of stars in the [C/Fe] — [Fe/H] plane is shown in Figure 3. 
The star-to-star scatter is large and a significant fraction of the model stars display a carbon 
enhancement in excess of [C/Fe] = 1. Since low-mass stars presumably were unable to 
form in gas with [C/H] < —3.5, hyper metal-poor stars are expected to be strongly carbon- 
enhanced, as observed. The absence of stars with iron abundances —5 < [Fe/H] < —4 is 
seen as a depression in the density function. The size of this gap depends primarily on the 
SFR, as well as on Knix, cf. Eq. (3), while the number of (carbon-enhanced) stars in the 
peak(s) below [Fe/H] = —5 also depends on the fraction of black hole-forming SNe. Note 
that, irrespectively of iron abundance, the majority of low-mass stars below [Fe/H] = —3 
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were formed during the initial phase of the second star formation epoch. 

In the present model, nearly 100% of the ISM was affected by ejecta from intermediate- 
mass stars at z — 5. However, the amount of carbon produced in this type of stars is 
not nearly enough to enrich the ISM to the level required by observations. The highest 
abundances of carbon can only be found near recent explosion sites of massive, black hole- 
forming SNc. These stars produce a lot of carbon and if all of it escapes the black hole, 
the ambient medium should indeed exhibit high C/H ratios. Nonetheless, for the generic 
set of yields (Meynet & Maeder 2002, Figure 3a), the model cannot explain the presence 
of the most carbon-rich stars. As suggested in Figure 3b, this discrepancy could largely 
be due to uncertainties in the calculations of the stellar yields. Using a set of models with 
somewhat different input physics, Meynet, Ekstrom, & Maeder (2006) found that the carbon 
yield increased a factor of 4 with respect to their generic calculations, at least for the most 
massive stars. It should also be noted that ID metal-poor model atmospheres overestimate 
carbon abundances derived from molecular lines with as much as 0.6 — 0.7 dex (see Collet, 
Asplund, & Trampedach 2005; Frebel et al. 2006). 

The scenario presented here predicts the possible existence of a small group of mega 
metal-poor stars ([Fe/H] < —6, see Figure 3). These stars would be enriched by a Pop III 
0-Ne-Mg SN (plus > 1 second generation black hole-forming SN). If so, as many as one out 
of five hyper metal-poor stars could be a mega metal-poor star. Their precise location in the 
[C/Fe] — [Fe/H] diagram is, however, uncertain due to the uncertain Fe yield of 0-Ne-Mg 
SNe. Not visible in Figure 3 is a second group of stars. These stars are formed out of gas 
only enriched by material from massive, second generation black hole-forming SNe. Although 
they are effectively iron-free, they may have a carbon abundance as high as [C/H] ~ —1.5, 
similar to the mega and hyper metal-poor stars. The predicted fraction of Fe-free stars in 
the Galactic halo is estimated to ~ 10""^. This fraction, however, depends sensitively on the 
SFR and the adopted primordial IMF (e.g., Nakamura & Umemura 2001; Abel, Bryan, & 
Norman 2002). 

4. Discussion and Conclusions 

An alternative source to the pre-enrichment of the Galactic halo may be the elusive pair- 
instability SNe, whose very massive progenitors (140 < M./M.Q < 260) have been suggested 
to dominate the first generation of stars (Bromm, Coppi, & Larson 1999; Abel et al. 2002; 
but see, e.g.. Silk & Langer 2006). If very massive stars indeed were formed in collapsing 
mini-halos surrounding the Galaxy, the reason why we have yet not found any trace of their 
nucleosynthetic signature in stars above [Fe/H] = —4 (see e.g., Tumlinson, Venkatesan, & 
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ShuU 2004) could thus be due to the fact that their ejecta were highly diluted and mixed 

with much less diluted ejecta of normal core-collapse SNe before locked up in low-mass stars. 
If so, the only observable signature of pair-instability SNe would be imprinted in the iron-peak 
of hyper metal-poor stars. It is, however, not clear whether their ejecta were able to mix 
with enough mass of interstellar gas (~ 10^ M.q) prior to the collapse of the Galactic halo. 

Apart from carbon, strong enhancements of nitrogen and oxygen are also measured in 
hyper metal-poor stars. For example, the [N/C] ratio in HE 0107—5240 and HE 1327—2326 
(subgiant solution) are —1.1 and +0.3, respectively (Frebel et al. 2006 and references 
therein). Non-rotating models of massive stars produce very little primary nitrogen (via 
the CNO-cycle) since freshly synthesized C and O never reach the H-burning shell (see, e.g., 
Umeda & Nomoto 2002; Meynet et al. 2006). In contrast, rotating (metal-poor) models 
are able to produce large amounts of nitrogen and indeed, the set of models (with different 
Vini and Z) by Meynet et al. (2006) encompasses the observed range in [N/C]. Rotating 
massive stars thus provide a plausible explanation to the enhancements of nitrogen in hyper 
metal-poor stars, in agreement with the above hypothesis. However, additional enrichment, 
such as mass-transfer from a possible binary companion, cannot be ruled out at this point. 
Regarding oxygen, a fraction of the hyper metal-poor stars in the simulations tend to be too 
oxygen-rich as compared to present observations. Since the oxygen yield depends strongly 
on stellar mass, the discrepancy can be reduced by assuming a lower limiting mass m„ for 
direct black hole-formation. However, a number of extremely metal-poor stars are display- 
ing oxygen abundances in excess of [0/H] = -2 (e.g., CS 29498-43, Aoki et al. 2004), in 
accordance with predictions. Thus, there may be hyper metal-poor stars with [0/H] > —2 
yet to be found. This merits further investigation. 

Observations of extremely and hyper metal-poor stars provide valuable information on 
the conditions at the time of formation of the Galaxy and it is imperative to, e.g., establish 
the metallicity distribution in these regimes. As discussed above, the lack of halo stars in 
the range —5 < [Fe/H] < —4 (Figure 2) suggests that the Galactic halo experienced a period 
of low star formation (star formation might have still occurred in e.g., the Galactic bulge) 
posterior to a prompt injection of metals by the first generation of stars. Moreover, the very 
strong enhancements of carbon, as well as of nitrogen and oxygen, in HE 0107—5240 and 
HE 1327—2326 are comparable to what is expected from subsequent enrichment of the ISM 
by rotating, black hole-forming SNe alone, without invoking additional sources (Figure 3). A 
significant fraction of the extremely metal-poor stars is also expected to be carbon-enhanced. 
Finally, if the inner and outer halo collapsed at different times, a correlation between iron 
abundance and orbital eccentricity of the hyper metal-poor stars may be expected to exist. 
If observed, such a correlation can put strong constraints on the formation of the Milky 
Way halo. Further modelling is, however, necessary in order to determine the level of pre- 
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enrichment and estimate effects of violent merging and density gradients in the ISM. 

I thank G. Meynet for discussions and for giving me a preprint of Meynet et al. (2006) in 
advance. I also thank A. Andersen, R. Collet, M. Davies, A. Prebel, B. Gustafsson, and E. 
Olsson, as well as the anonymous referee, for valuable comments and suggestions. 

REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2002, Sci, 295, 93 

Aoki, W., Norris, J. E., Ryan, S. G., Beers, T. C., & Ando, H. 2002, ApJ, 576, L141 

Aoki, W., et al. 2004, ApJ, 608, 971 

Barklem, R S., et al. 2005, A&A, 439, 129 

Bateman, N. P. & Larson, R. B. 1993, ApJ, 407, 634 

Beers, T. C. 1999, Ap&SS, 265, 105 

Beers, T. C. & Christlieb, N. 2005, ARA&A, 43, 531 

Bromm, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 

Bromm, V. & Loeb, A. 2003, Nat, 425, 812 

Christlieb, N., et al. 2002, Nat, 419, 904 

Collet, R., Asplund, M., & Trampedach, R. 2005, in lAU Symp. 228, From Lithium to 
Uranium: Elemental Tracers of Early Cosmic Evolution, ed. V. Hill, P. Prangois, & 
P. Primas (Cambridge: CUP), 247 

Daigne, P., Olive, K. A., Vangioni-Plam, E., Silk, J., & Audouze, J. 2004, ApJ, 617, 693 

Prebel, A., et al. 2005, Nat, 434, 871 

Prebel, A., Christlieb, N., Norris, J. E., Aoki, W., & Asplund, M. 2006, ApJ, in press 
(astro-ph/0512543) 

Garcia-Percz, A. E., Asplund, M., Primas, P., Nissen, P. E., & Gustafsson, B. 2006, A&A, 
in press (astro-ph/0512290) 

Iwamoto, N., Umeda, H., Tominaga, N., Nomoto, K., & Maeda, K. 2005, Sci, 309, 451 



-9- 

Karlsson, T. 2005, A&A, 439, 93 

Karlsson, T. & Gustafsson, B. 2005, A&A, 436, 879 

Kroupa, P. 2001, MNRAS, 322, 231 

Maylc, R. & Wilson, J. R. 1988, ApJ, 334, 909 

Meynet, G., Ekstrom, S., & Maeder, A. 2006, A&A, in press (astro-ph/05 10560) 

Meynet, G. & Maeder, A. 2002, A&A, 390, 561 

Nakamura, F. & Umemura, M. 2001, ApJ, 548, 19 

Plez, B. & Cohen, J. G. 2005, A&A, 434, 1117 

Reed, D. S., et al. 2005, MNRAS, 363, 393 

Ryan, S. G., Norris, J. E., & Beers, T. C. 1996, ApJ, 471, 254 

Samland, M., Hensler, G., & Theis, C. 1997, ApJ, 476, 544 

Silk, J. & Langer, M. 2006, MNRAS, submitted 

Sommer-Larsen, J., Gotz, M., & Portinari, L. 2003, ApJ, 596, 47 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Spite, M., et al. 2005, A&A, 430, 655 

Springel, V. & Hernquist, L. 2003, MNRAS, 339, 312 

Tumlinson, J., Venkatesan, A., & ShuU, J. M. 2004, ApJ, 612, 602 

Umeda, H. & Nomoto, K. 2002, ApJ, 565, 385 

— . 2003, Nat, 422, 871 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 10 - 



Time (Myr) 




25 20 15 10 

Redshift, z 



Fig. 1. — The adopted SFR in the Galactic halo (black line) as given by the fraction of 
massive stars exploding as SNe. The shaded area indicates the period of low SFR. 
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Fig. 2. — The predicted halo metalhcity distribution. The dashed hne denotes the distribu- 
tion of stars neglecting the suppression of low-mass star formation in carbon-deficient gas. 
The full line denotes the fraction of stars with carbon abundance [C/H] > — 3.5. In both 
cases, a distinct population of stars appears in the range —5.5 < [Fe/H] < —5, in agreement 
with observations (symbols). 
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Fig. 3. — The predicted distribution of stars in the [C/Fe] — [Fe/H] plane. Upper panel: The 
thick, dashed hne indicates a carbon abundance of [C/H] = —3.5. The shaded area below 
this limit should contain no low- mass stars (cf. Figure 2). The symbols denote observations 
of stars in the Milky way's halo. The circles are unevolved dwarf/subgiant stars (Barklem 
et al. 2005), the diamonds are stars on the lower red giant branch (Spite et al. 2005), the 
square is G 77—61 (Plez & Cohen 2005), the upward facing triangle is CS 29498—43 (Aoki 
et al. 2004), the downward facing triangle is CS 22957-027 (Aoki et al. 2002), and finally, 
the pentagon and hexagon are HE 0107-5240 (Christlieb et al. 2002) and HE 1327-2326 
(Frebel et al. 2005), respectively. G 77—61 and CS 22957—027 are known to be members of 
binary systems and may have been born with a lower surface carbon abundance as indicated 
by the arrows. The predicted group of mega metal-poor stars is encircled. Bottom panel: 
same as above but with a carbon yield increased by a factor of 4 for stars in the mass range 
30 < m/A4Q < 60. Note that observed abundance ratios are not corrected for 3D effects. 



